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^ ! Abstract 

We present quasi-elastic neutrino-nucleus cross sections in the energy range from 150 MeV up 
CN ' to 5 GeV for the target nuclei 12 C and 56 Fe. A relativistic description of the nuclear dynamics and 

> ; 

the neutrino-nucleus coupling is adopted. For the treatment of final-state interactions (FSI) we 



o 
in 
o 
in 
o 



rely on two frameworks succesfully applied to exclusive electron-nucleus scattering: a relativistic 
optical potential and a relativistic multiple-scattering Glauber approximation. At lower energies, 
the optical-potential approach is considered to be the optimum choice, whereas at high energies 
a Glauber approach is more natural. Comparing the results of both calculations, it is found that 
^ \ the Glauber approach yields valid results down to the remarkably small nucleon kinetic energies 



of 200 MeV. We argue that the nuclear transparencies extracted from A(e, e'p) measurements can 



^ . be used to obtain realistic estimates of the effect of FSI mechanisms on quasi-elastic neutrino- 

h : 

nucleus cross sections. We present two independent relativistic plane-wave impulse approximation 
(RPWIA) calculations of quasi-elastic neutrino- nucleus cross sections. They agree at the percent 
level, showing the reliability of the numerical techniques adopted and providing benchmark RPWIA 
results. 

PACS numbers: 25.30.Pt; 13.15.+g; 24.10.Jv 
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I. INTRODUCTION 



Neutrino interactions offer unique opportunities for exploring fundamental questions in 
different domains of physics. The mass of the neutrino remains one of the greatest puzzles 
in elementary particle physics. In recent years, a number of positive neutrino oscillation 
signals made the claims of non-zero neutrino masses irrefutable and boosted the interest 
in this issue. Several experiments are running or proposed in order to address intriguing 
questions in current neutrino physics Q|: What does the neutrino mass hierarchy look like, 
and what are the values of the oscillation parameters lj? What is the role of vacuum and 
matter-enhanced oscillations ? Are neutrinos representatives of CP-violation in the leptonic 
sector ? Is the neutrino a Dirac or a Majorana particle ? Does it have a magnetic moment 

? s. 

The interest in neutrinos goes beyond the study of the particle's intrinsic properties, and 
extends to a variety of topics in astro-, nuclear and hadronic physics. Typical astrophysical 
examples include the understanding of the energy production in our sun, neutrino nucleosyn- 
thesis and the synthesis of heavy elements during the r-process, the influence of neutrinos on 
;he dynamics of a core-collapse supernova explosion and the cooling of a proto-neutronstar 
. In many astrophysical situations the neutrinos serve as messengers probing the interior 
of dense and opaque objects that otherwise remain inaccessible. The influence of neutrinos 
even extends to cosmological questions such as the role of neutrinos in the matter-antimatter 
asymmetry in the universe. In hadronic and nuclear physics, neutrino scattering can shed 
light on a lot of issues, including investigations of electroweak form factors, the study of the 
strange quark content of the nucleon and ^-induced pion production 
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Despite the richness of phenomena they are involved in, neutrinos remain elusive par- 
ticles, only weakly interacting and eager to escape detectors on the watch. The presence 
of neutrinos, being chargeless particles, can only be inferred by detecting the secondary 
particles they create when colliding and interacting with matter. Nuclei are often used as 
neutrino detectors, providing relatively large cross sections that offer a broad variety of in- 
formation. As a consequence, a reliable interpretation of data involving neutrinos heavily 
counts on a detailed knowledge of the magnitude of neutrino-nucleus interactions under var- 
ious circumstances. A precise knowledge of the energy and mass number dependence of the 
neutrino-nucleus cross section is essential to current and future measurements. The energies 
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that neutrinos can transfer to nuclei depend on their origin. The 'low' energy regime extends 
to a few tens of MeV and relates to reactor, solar and supernova neutrinos. Atmospheric 
and accelerator neutrinos can carry energies from a few hundred MeV to several GeV's. 

At intermediate energies (here defined as energies beyond the nuclear resonance region), 
neutrino-nucleus interactions have been studied within several approaches, investigating a 
variety of effects. The relativistic Fermi gas (RFG) model was employed in Refs. ja, E| to 
study the possibility of measuring strange-quark contributions to the nucleon form factors. 
The RFG takes into account the Fermi motion of the nucleons inside the nucleus, Pauli block- 
ing and —tic Katies, bu t ne gl ect s severa, other effect, Re, Q used a pla_ 
impulse approximation description of the nuclear system to estimate polarization asymmetry 
effects in neutrino-induced nucleon knockout. Relativistic nuclear effects were included in 
the calculations of Refs. Q, Q, Q, , 
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171 ]. using a relativistic shell model approach 



for the study of neutral-current and/or charged- current neutrino-nucleus scattering. In par- 
ticular, in Refs. Q] results in the relativistic plane-wave impulse approximation 
(RPWIA) were compared to RFG calculations. It is shown that binding-energy effects tend 
to vanish as the energy increases. Going one step further in the complexity of the model 
calculation, the implementation of the final-state interactions (FSI) of the ejected nucleon 
has been achieved in different manners. In Ref. a phenomenological convolution model 
was applied to the RFG, showing that nucleon re-scattering can produce a reduction of the 
quasi-elastic cross section as large as 15% at incoming neutrino energies of about 1 GeV. 
A description of FSI mechanisms through the inclusion of relativistic optical potentials is 
presented in Refs. More specifically, Ref. |1J] studies the uncertainties 
derived from the use of different prescriptions for the potentials. A reduction of the cross 
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section of at least 14% is found at incoming neutrino energies of 1 GeV. In Refs. 
important FSI effects arise from the use of relativistic optical potentials within a relativistic 
Green's function approach. Apart from relativistic dynamics and FSI, other effects may 
have an impact on neutrino-nucleus reactions. In Refs. jlfll . I^fll ] the influence of relativistic 
nuclear structure effects, delta- and pion degrees-of-freedom, and RPA-type correlations on 



neutrino-scattering cross sections was examined. Ref. [2l| includes long-range correlations, 
FSI and Coulomb corrections in 12 C(^, fi~) 12 C* calculations. An alternative method was 
proposed in Ref. J%J, where it was shown that a superscaling analysis of few-GeV inclusive 
electron scattering data allows one to predict charged- current neutrino cross sections in the 
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nuclear resonance region, thereby effectively including delta isobar degrees-of-freedom. 

In this paper we compute the single-nucleon knockout (often referred to as quasi-elastic 
(QE)) contribution to the inclusive neutrino- nucleus cross sections, for energies and nu- 

I. We consider 



clei relevant to proposals like Minerva Q], Miniboone [f| and Finesse 
that the large variety of relevant neutrino energies and the tendency to study neutrino- 
nucleus interactions at increasing energies, necessitate the use of relativity. We employ 
two relativistic models for describing neutrino-nucleus scattering within the impulse ap- 
proximation: the relativistic distorted-wave impulse approximation (RDWIA) developed by 
the Madrid-Seville group, and the relativistic multiple scattering Glauber approximation 
(RMSGA) developed by the Ghent group. Initially designed for the description of exclu- 
sive electron-nucleus scatteri ng p rocesses, both models have been succesfully tested against 



A(e,e'p) data [2J, |25j, |2y, |27|, |28|, |29|, |30| • In addition, the nuclear transparencies predicted 
by these models have proven to be mutually consistent in the intermediate kinematic re gim e 



^gime 



between 0.5 and 1 GeV nucleon kinetic energies where both of them are deemed reliable 
The RDWIA model used here has already been employed in several neutrino-nucleus cal- 
culations Q]. To our knowledge, this paper is the first report of a relativistic 
Glauber-inspired approach to neutrino-nucleus reactions. The aim of this work is twofold. 
First, the relativistic models available to date predict different results in the limit of van- 
ishing FSI, motivating a 'new round' of calculations. We investigate the plane-wave limit of 
the RDWIA and RMSGA approximations, aiming at providing benchmark RPWIA results. 
Second, we compute the effects of FSI within our models, paying special attention to the 
comparison between RDWIA and RMSGA results. It is well known that the inclusion of FSI 
within inclusive calculations requires a considerable computing effort. We propose a way to 
estimate FSI effects for the QE contribution to the inclusive neutrino-nucleus cross section 
using benchmark RPWIA results and transparency data obtained from A(e, e'p) experi- 
ments. For the time being, the effects of many-body currents, nucleon-nucleon correlations, 
and contributions beyond quasi-elastic scattering processes as multi-nucleon processes and 
pion production are neglected. 

The outline of this paper is as follows. In Sec. |H]we present the RMSGA and RDWIA for- 
malisms for the description of the neutral- and charged-current neutrino-nucleus scattering 
processes. Results of the numerical calculations are shown in Sec. IIIII Sec. IIVI summarizes 
our findings. 



4 




q"=(a>,q) 



^A-l = ^A-l*^A-li 



K$=(E f ,k f ) k 
A / / 



FIG. 1: Kinematics for the quasi-elastic neutrino- nucleus scattering process 



II. FORMALISM 

We derive expressions for neutrino and antineutrino neutral-current (NC) reactions from 
nuclei which result in one emitted nucleon 



Here, / labels the flavor of the lepton, and A represents a nucleus with mass number A. 
The connection between electromagnetic and weak interactions makes that the analytical 
derivations follow the same lines as those used in electron-nucleus scattering. The main 
differences between neutrino and electron interactions stem from the intrinsic polarization 
of the neutrino due to the parity-violating character of the weak interaction. Moreover in 
weak interactions the focus is on inclusive processes, whereas exclusive processes play a 
predominant role in current subatomic research with electrons. 

We describe these processes at lowest order in the electroweak interaction, i.e. considering 
the exchange of one charged vector boson. Fig. ^defines our conventions for the kinematical 
variables. The four-momenta of the incident neutrino and scattered lepton are labeled 
and K' 11 . Further, Kj^_ x and represent the four-momenta of the target nucleus, the 




(1) 



We also consider their charged-current (CC) counterparts 




(2) 
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residual nucleus and the ejected nucleon. The xyz coordinate system is chosen such that 
the z— axis lies along the momentum transfer q, the y— axis along k x k' and the x— axis lies 
in the scattering plane. The hadron reaction plane is then defined by kf and q. We adopt 
the standard convention Q 2 = —q^ for the four-momentum transfer. 



A. Quasi-elastic neutrino-nucleus cross section 



In the laboratory frame, the exclusive differential cross section for the processes specified 
in Eqs. (JTJ) and © can be written as 



da = -£|M 



M x M A -xM_ 



fi\ 



Nj3 



if 



e' E A -i E 



f 



d A k A -id 6 k'd A k s 



(3) 



where Ylij indicates sum and/or average over initial and final spins. Dealing with neutrinos, 
the relative initial velocity j3 can trivially be put to 1. The factor ^ stems from the 
normalization of the outgoing lepton spinor and becomes 1 for NC reactions. Integrating 
over the unobserved momentum of the recoiling nucleus fc^-i, as well as over \kf\, results 
in the following fivefold differential cross section for the A(u, v'N), A(V,V'N), A(u, IN) and 
A(u,lN) reactions 



d 5 a 



de'd^d^f (2ir) 5 M A e 



M,M N M Afk%f _^ Mf 
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(4) 



if 



where Qi and Qf define the scattering direction of the outgoing lepton and the outgoing 
nucleon. The recoil factor f rec is given by 



fr 



E A - 



Ma 



1 + 



E 



■[1- 



A-l 



(5) 



The squared invariant matrix element Mjj can be written as 



if 



Ml 



Q 2 



(6) 



Here, Mb represents the mass of the Z-boson for NC reactions and that of the IV-boson for 
CC processes. Gp is the Fermi constant. For CC reactions the latter has to be multiplied 
with a factor cos^, with 9 C the Cabbibo angle, determining the mixing of the strong down 
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and strange quarks into the weak d-quark. In the above expression the lepton tensor is 
defined as 

la/3 = ^[Ullail ~ l5)ui] j [u u lp(l ~ l5)Uu], (7) 
s,s' 

with s and s' the initial and final lepton spins. The hadron tensor is given by 

W af > = ^(A^J^A^) = EUlV"), (8) 

if if 



with the boson propagator 



q n qV 



l B 

The quantity (J~ a ) in Eq. (JBJ) can be written as 



A"" = f v - ^r- (9) 



(J a ) = \{A- 1)(J R M R ), K f (E f , k f )m s A afl J» A(0 + ,g.s.)j, (10) 

with J M the weak current operator, \A(0 + , g.s.)) the ground state of the target even-even 
nucleus and \ (A — 1)(JrMr)) the state in which the residual nucleus is left. At the energies 
considered here, Eq. (JHJ can approximately be written as A M1/ m g^, and the quantity 
{J a ) ~ {J a )i establishing the connection between the four- vector J a and the nuclear current 
operator. In the extreme relativistic limit, the contraction of the lepton tensor l a p with the 
nuclear one W a/3 in Eq. © can be cast in the form : 



d 5 a MnMa-i ,_i z,w± 

de'd 2 n i( Pn f ~ (2ir) S M A fh ec°M 



with a M defined by 



for NC reactions and 



x [vlRl + vtRt + vttRtt COS 2(f) 

+vtlRtl cos0 + h{v' T R' T + v' TL R' TL cos0)] , (11) 



^ ,' GFCQs(V2KMf V 



V2-k{Q 2 + M : 



° M = V 1 -^U-(Q 2 + M^) ) ' (13) 

for CC reactions. In these equations, Q\ is the angle between the direction of the incident 
and the scattered lepton's momentum and <p the azimuthal angle of the reaction plane (see 
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Fig. H}. In Eq. h = — 1 (h = +1) corresponds to the helicity of the incident neutrino 
(antineutrino). For NC reactions, the lepton kinematics is contained in the kinematic factors 



v L = 1, (14) 
v T = tan — + —r^, (15) 



J TL 



VTT = ~W\ (16) 

VTL = ~71V * W' (17) 



taa|jtaa 2 | + ^, (18) 



2 V 2 l<?1 



^tan|. (19) 



The corresponding response functions read 



Rt = \(J + (q))\ 2 + \(J~{q))\\ 
R TT cos2<p = 2$t{(J + (q))*(J-(q))}, 

R TL cos4> = -2^{\(J°(q))-^(J°(q)) 
I L 1 91 

R't = \(J + m\ 2 -\{J-{q))\\ 



R' TL cos ( 



-23? 



UJ 



(J\q))-^(J z {q)) 



[(J + (q))-{J-(q))Y 



where {J(q)) is expanded in terms of unit spherical vectors e r 



e 



e±i = =f-j=(e x ± ie y ). 



(20) 

(21) 
(22) 
(23) 
(24) 
(25) 



(26) 



For CC reactions, the mass of the outgoing lepton has to be taken into account. This 
results in the following substitutions (see also for instance j^|) 
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VTT 



vtl 



J TL 



c'2 



COS 01 



EE' 
W 



Ml 

all 



sin 2 0i, 



EE' 

sin 0i 



Mt 



sin 2 0i 



e + e') 



Furthermore 

Rtl cos< 

and 



23? 



1*1 

sin 0[ 

'^7T 



i - 



COS 01 



Ml 



Ml 



\q\ 



[(J + {q))-(J-{q))X 



v L R L = v° L R L + vlRl + v L z R° L *, 



with 



and 



Ri = \(j (q))\ 2 , Ri = \{J z m\- 



-2St{(J°(q))(J z Wr} 



M i a 



.02 



1 + A 



M? 



-12 



COS 01 — 



2ee' 



1 - 



M, 



Ml 

c '2 



COS 01 + — f 



sin 2 ^ 



M 2 



(27) 
(28) 
(29) 

(30) 
(31) 



(32) 

(33) 
(34) 

(35) 
(36) 
(37) 



The expressions for R T , Rtt, R't an d R'tl remain unaltered. 

So far, a precise knowledge of the kinematic variables at the lepton vertex was assumed. 
In practice, this information is not attainable in typical neutrino scattering experiments. 
Indeed, in NC reactions, the scattered lepton is chargeless and remains undetected. In 
CC processes, on the other hand, detection of the final lepton is possible and its energy 
and momentum could in principle be measured. However, due to limited control on the 
incoming neutrino energies, the energy-momentum balance at the lepton vertex cannot be 
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precisely determined. In order to get the QE neutrino-nucleus cross section, we integrate 
over the phase space of the scattered lepton (d 2 Qi) and the outgoing nucleon (d 2 flf(9f, </>)). 
For the latter, integration over the azimuthal angle <fi yields a factor 2tt, whilst only the 
(^-independent terms survive due to symmetry properties. This yields 



In practice, we compute the response functions for all single-particle levels in the target 
nucleus, and obtain dajde' by summing over all these. 

B. Nuclear current 

Obviously, the determination of the response functions requires knowledge of the nuclear 
current matrix elements tjlOj) . We describe the neutrino- nucleus nucleon- knockout reaction 
within the impulse approximation, assuming that the incident neutrino interacts with only 
one nucleon, which is subsequently emitted. The nuclear current is written as a sum of single- 
nucleon currents. The wave functions for the target and the residual nuclei are described in 
terms of an independent-particle model. Then, the transition matrix elements can be cast 
in the following form : 



where 4>b and <\>p are relativistic bound-state and scattering wave functions. Further, J M is 
the relativistic one-body current operator modeling the coupling between the virtual Z° or 
W ± boson and a bound nucleon. The relativistic bound-state wave functions are obtained 
within the Hartree approximation to the <j-oj model j^. The quantum- field theoretical 
problem is solved in the standard mean-field approximation replacing the meson field op- 
erators by their expectation values. The resulting eigenvalue equations of the relativistic 
mean-field theory can be solved exactly. The corresponding bound-state wave functions 4>b 
are four-spinors and can formally be written as 




\d6i \ sin 9 f d9 fk f f re lcTM[vLRL + v t Rt + hv' T R' T } . 



(38) 




(39) 




(40) 



with 3^ s m B (^r, g) the familiar spin spherical harmonics. 
We use a relativistic one-body vertex function of the form: 
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J» = F 1 (Q' 2 )r + + G A (Q 2 )Yl 5 + ^^(QVts, ( 41 ) 

with k the anomalous magnetic moment. The weak vector form factors F\ and F2 can be 
related to the corresponding electromagnetic ones for protons (Ff^ 1 ) and neutrons (FfJ^) 
by the conserved vector current (CVC) hypothesis. For proton knockout they are given by 

(i - 2 sin 2 6 W ) F? p M - \F™ for NC reactions, 

{F? p M - FF? 1 ) for CC reactions, 

with 9w the Weinberg angle defined by sin 2 9\y = 0.2224. For neutron knockout the weak 
vector form factors result from the exchange of the subindexes p and n in Eq. (|42j). A 
standard dipole parametrization is adopted for the vector form factors. 
The axial form factor for proton knockout is expressed as 

-M G for NC reactions, 

(43) 

qaG for CC reactions 

where g A =1.262, an d G = (1 + Q 2 /Af*)~ 2 with M = 1.032 GeV. For neutron knockout a 
minus sign must be added to Eq. 

The Goldberger-Treiman relation allows one to write the pseudoscalar form factor as 

Gp(Q 2 ) = ^^G A (Q% (44) 

where m, denotes the pion mass. The contribution of this form factor, being proportional 
to the mass of the scattered lepton, vanishes for NC reactions. 



C. Final-state interactions in relativistic models 



We now turn to the question of computing a relativistic scattering wave function for 
the outgoing nucleon. Including nucleon-nucleus FSI is a long-standing issue in theoretical 
A(e, e'p) investigations. For kinetic energies up to around 1 GeV, most calculations have 
traditionally been performed within a so-called distorted-wave impulse approximation model 
(DWIA), where the final nucleon scattering state is computed with the aid of proton- nucleus 
optical potentials. For proton kinetic energies above 1 GeV parameterizations of these 
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potentials within the context of Dirac phenomenology are not readily at hand. Furthermore, 
beyond this energy the use of optical potentials for modeling FSI processes does not seem 
very natural in view of the highly inelastic and diffractive properties of the underlying 
nucleon-nucleon scattering process. In this energy regime, the Glauber model, which is a 
multiple-scattering extension of the eikonal approximation, offers a valid and economical 
al—e for descriW FSI Q. In a Glanner t _ k , the FSI effect s are computed 
directly from the elementary nucleon-nucleon scattering data. Below, we give a brief outline 



of the main features of both mode 
Within the RDWIA framework 



s. 

24Ll25l.l2fil.l27l]. (pp in Eq. (jHHJ) is a scattering solution to a 



Dirac-like equation, which includes scalar and vector complex optical potentials obtained by 
fitting elastic pA scattering data. The real part of these potentials describes the rescatterings 
of the ejected nucleon. The imaginary part accounts for the absorption into unobserved 
channels. The scattering wave function, expressed in terms of a partial-wave expansion in 
configuration space, reads 



J K/j,m 

xYr(n kf )K(?), (45) 



where \&£(r) are four-spinors of the same form as in Eq. (J4~U|) . The phase-shifts 5* and 
radial functions are complex because of the complex potentials. The outgoing nucleon spin 
is denoted as sj. In this work we use the relativistic global optical potential corresponding 
to the energy and target mass-dependent parametrization (EDAD1) of Ref. 3sj |. 

The Glauber approach relies on the eikonal and the frozen-spectators approximation. 
It allows to formulate a full-fledged multiple-scattering theory for the emission of a "fast" 
nucleon from a composite system consisting of A— 1 temporarily "frozen" nucleons. Ref. 
provides a detailed outline of a relativistic and unfactorized formulation of Glauber multiple 
scattering theory. In this approach, coined RMSGA, the scattering wave function in the 
matrix element of Eq. (pHUj) adopts the form 

<f> F (^)=G(b,z)<f> kf>Sf (f) (46) 

where <pk f , Sf is a relativistic plane wave. The impact of the FSI mechanisms on the scattering 
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wave function is contained in the scalar Dirac-Glauber phase Q(b, z) 
G(b,z)= J] I- j dr>\<j> a {r')\%z' -z)T(jb' -b) 



(47) 



where the product over a(n,n,m) extends over all occupied single-particle states in the 
target nucleus, excluding the one from which the nucleon is ejected. The profile function for 
NN scattering is defined in the standard manner 



r(6) 



a 



tot I 
NN ' 



ie NN ) 



■ CX P( 



^Pnn ^Pnn 



(48) 



The parameters <r^v, /3nn and e^N depend on the ejectile energy. Values for the parameters 
fitted to the pp and pn data can be found in Ref. 36j |. The neutron-neutron scattering 
parameters are assumed identical to the proton-proton ones. 

As the integrations in Eq. (|38|) would require an enormous numerical effort, we introduce 
an additional averaging over the positions of the spectator nucleons. This procedure amounts 
to replacing in Eq. (|4Tj) the characteristic spatial distributions of each of the spectator 
nucleons by an average density distribution for the target nucleus 



ItNN. 



X 



^Pnn 
b'db'T B (b',z)exp 



A2 



(b ~ V) 

2(3nn 



2tt 



%i exp 



-2W 



NN 



■sm 



A-l 



(49) 



The function Tsib', z) which was introduced in the above expressions is known as the "thick- 
ness function" and reads 

i r + °° 

T B (b', z) = - dz'9{z' - z)p B (r'(b', z')), (50) 

A J-oo 

where the relativistic radial baryon density p_e(r) is defined in the standard fashion 
p B (r) = £ / dadn{4> a {r,G))\Ur,a)) 

a J 

(2j + l 



E 



Airr 2 



\G n Jr) 



\Fn*(r)\' 



(51) 



and the sum over nn extends over all occupied states. For exclusive A(e,e',p) processes, 
where the quantum numbers of the residual nucleus are well defined, the thickness-function 
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approximation of Eq. (j 
expression of Eq. (JITj) 



provided results which approach the exact ones obtained with the 
|. Here, we only deal with inclusive cross sections, obtained after 
incoherently summing over nucleon emission from all possible single-particle shells. Under 
these circumstances, one can expect that the thickness-function approximation becomes an 
even better one. 

The RDWIA and RMSGA codes were developed independently and adopt distinctive 
numerical techniques to compute the scattering wave functions and the corresponding ma- 
trix elements of Eq. (|3*^j) . The RDWIA code employs a partial- wave expansion to solve the 
Dirac equation for the ejectile. The cylindrical symmetry of the Glauber phase of Eq. (|47|) 
prohibits any meaningful use of this technique in the RMSGA calculations. Instead, the 
multi-dimensional integrals are computed numerically. In the limit of vanishing FSI mech- 
anisms, i.e. within RPWIA, though, both models should yield identical results. In the 
Glauber picture this limit is reached by putting the Glauber phase of Eq. (j4Tlj) equal to 
unity. In the RDWIA picture, the effect of FSI can be made vanishing by nullifying the 
optical potentials. Then, the computed partial waves sum to a relativistic plane wave. Con- 
vergence of the partial wave expansion was tested against the analytical plane- wave result. 

The models described above were initially developed for the description of exclusive 



A(e, e'p) processes, for which an exce 



data has been achieved (24 , 



25 



26 



27 



lent agreement between theoretical calculations and 
j. It is clear that inclusive neutrino scattering 



cross sections include contributions which fall beyond the scope of the RDWIA and RMSGA 
models. Both the RDWIA and RMSGA are confined to those processes where the scattering 
of a neutrino from a nucleus causes a single nucleon to escape, thereby exciting the residual 
nucleus 111 cl StcltG tit missing energies below 80 MeV and a predominant single-hole nuclear 
structure with respect to the ground state of the target nucleus. We refer to such processes 
as "elastic" ones and wish to stress that they include proton and neutron knockout from 
the deepest lying Is up to the Fermi level. Inelastic single-nucleon knockout channels pop- 
ulating more complex states in the residual A-l nucleus are excluded from our calculations. 
So are multi-nucleon knockout channels and channels involving a pion. In that sense, the 
RDWIA and RMSGA predictions for the inclusive neutrino-nucleus cross sections should be 
interpreted as a lower limit of the single-nucleon knockout contribution. 
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III. RESULTS 



We present results for QE neutrino scattering from 12 C and 56 Fe, which are nuclei well 
suited for neutrino detection. The calculations span incident neutrino energies from 150 MeV 
up to 5 GeV. From about 200 MeV to 1 GeV, the quasi-elastic nucleon knockout is expected 
to be the dominant contribution to the neutrino-nucleus cross section. At higher energies, 
the relative contribution of the inelastic channels, mainly those involving an intermediate 
delta resonance a nd p ion production, is expected to become increasingly dominant in the 
inclusive process 

HQ- Ref. Q indicates that in the neutrino energy range from 0.7 to 
5 GeV reaction channels involving a pion contribute for 15% to the total cross section. 

In order to make the comparisons between the RDWIA and RMSGA calculations as 
meaningful as possible, all the ingredients in the calculations not related to FSI, as those 
concerning the implementation of relativistic dynamics and nuclear recoil effects, are kept 
identical. In particular, both pictures adopt the Wl parametrization for the different 
field strenghts in determining the bound-state wave functions. Accordingly, the RDWIA 
and RMSGA only differ in their assumptions regarding the treatment of FSI. 

It speaks for itself that before embarking on the study of effects like FSI, pion production, 
the role of the delta in the medium, multi-nucleon knock-out, the strangeness content of the 
nucleon it is absolutely essential to possess reliable baseline RPWIA cross sections with 
a numerical accuracy of a few percent. To this purpose, before turning to the study of the 
role of FSI mechanisms, we first investigate the RPWIA limit of the RMSGA and RDWIA 
models. These predictions will be compared and confronted with other RPWIA results 



which made their way to literature recently 
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A. RPWIA 

Fig.Elshows the results of various RPWIA calculations for 12 C(z/, z/) at 150, 500 and 1000 
MeV. We observe that the plane-wave limits of our RMSGA and RDWIA formalisms are in 
excellent agreement. The remaining differences, smaller than 2-3%, can be attributed to the 
distinctive numerical techniques. This comparison lends us confidence about the consistency 
of the two types of calculations and the reliability of the adopted numerical techniques. 

The fact that our models provide almost identical RPWIA results may seem trivial. As 
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FIG. 2: Neutral current 12 C(v, v') cross sections as a function of the outgoing nucleon kinetic energy 
Tjv at different incoming neutrino energies. The solid (dashed) lines represent the RPWIA results 



Q. 



of the Ghent (Madrid) group. The short-dot-dashed lines show the RPWIA results of Ref. 
and the long-dotted lines those of Ref. 2). The short-dotted (long-dot-dashed) line shows the 
predictions of the RFG model of Ref. 14] (Ref. ^|) with a binding-energy correction of 27 MeV. 



can be appreciated via Fig. |2J however, our RPWIA predictions disagree with the ones of 
Refs. [13 and 1?J. Although the RPWIA calculations of Refs. jl^ and [itJ are mutually 
consistent at e = 500 MeV, this is no longer the case at e = 1 GeV. In the search for the origin 
of the discrepancies between our and other RPWIA calculations, differences in the nuclear 



1(3 



current can be ruled out. The current operator of Eq. (|41|) used along this work is formally 
identical to the one mentioned in Refs. [151 and Il7|, and the same holds for the form- factor 

n n 

parameterization. Only the bound-state wave functions used in Refs. |15J and [17J differ 
from ours. We have performed cross section calculations with various parameterizations for 
the bound state wave functions, and found almost negligible differences. 

The role of the various terms F\, Fo and Ga in Eq. (HT1) in the NC differential cross 

. The results were shown for proton knockout from the 
lp 3 / 2 orbital of 12 C, at incident neutrino energies of 150, 500 and 1000 MeV. In Fig. El we 
analyse the contribution of the F\, F2 and Ga form factors in our cross sections under the 
same circumstances. First, let us observe that the calculations performed by setting F\ = 
(long-dot-dashed lines) almost reproduce the full cross sections (solid lines). This shows, in 
agreement with the outcomes of Ref. [l^ |. that the contribution of the Dirac form factor 
is very small. The cross section can then be very well approximated as a sum of three 
terms: one proportional to (Ga) 2 , a second proportional to (-F2) 2 , and a third proportional 
to the interference of Ga and F2 contributions. The term proportional to (Ga) 2 (dashed 
lines) is very similar to the corresponding one in Fig. 11 of Ref. ^J]. The same holds 
for the cross sections obtained by nullifying Ga (short-dot-dashed lines), whose behaviour 
is almost entirely given by F| term. Thus, the fact that our curves neglecting the Fi 
contribution differ from those in Ref. 3] can mainly be attributed to the GaFi interference 
term. Furthermore, the results seem to differ in the sign of this term. As a matter of fact, 
changing the sign of the term in our calculations, the cross sections follow closely the 



ones of Ref. Il7|. At 500 and 1000 MeV the differential cross sections of Ref. 



171 ] display 



some oscillations as a function of T^. As can be appreciated from Fig. El and H3 we find 
no indications for these oscillations. Recently, the authors of Ref. 17j have extended their 
work to calculate CC neutrino cross sections j4l|. We remark that for this type of neutrino 
reactions the magnitude of the cross sections is in agreement with the RPWIA limit of the 
models presented here. 

It is well known that binding-energy effects tend to vanish with increasing energies. 
Accordingly, a description of the ^-nucleus scattering process in terms of a RFG model is 
expected to approach the RPWIA predictions at high incoming, neutrino energies. This is 



observed in Fig. El when comparing the RFG results of Refs. [a |l2|, |lj] with our RPWIA 



predictions. At e — 150 MeV, our RPWIA cross sections are approximately 15% smaller 
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FIG. 3: Effect of the different form factors on the neutral current 12 C(u, u f ) cross sections as a 
function of T/v at different incoming neutrino energies. The solid lines represent the full RPWIA 
results. The long-dot-dashed (short-dot-dashed) lines show the results with F\ = {Ga = 0). The 
dashed lines show the cross section when only Ga is considered. 



than the RFG ones. The RFG result from Refs. |12t Il4j closely follows our RPWIA results 
at 500 MeV, the agreement at 1 GeV being remarkabl y g ood. The observed similarity 
between the independent RFG predictions of Refs. js, 12, 1A\ and our RPWIA results lends 
us additional confidence that the RPWIA results presented here can serve as benchmark 
calculations. 
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B. The effect of FSI : RMSGA and RDWIA approaches 



Let us now turn our attention to the effect of FSI. NC //-nucleus cross sections obtained 
within RDWIA and RMSGA are displayed in Fig. 0] The calculations correspond to 12 C 
and 56 Fe targets, and incoming energies of 500, 1000, and 5000 MeV. Focusing on the results 
of the RDWIA model, the inclusion of the complex optical potential reduces the RPWIA 
results by nearly 40 — 50% for 12 C. As expected, the global effect of FSI increases with 
growing atomic number, and reductions of over 60% are obtained for 56 Fe. The presence 
of the imaginary term in the optical potential is likely to lead to an underestimation of 
the single-nucleon knockout contribution to the inclusive cross section. Indeed, in inclusive 
measurements all possible final channels are included, whilst the RDWIA and RMSGA 
calculations are confined to "elastic" single-nucleon knockout. 

Traditionally, Glauber-inspired models have been esteemed to provide reliable results at 
high energies, as they rely on the eikonal approximation. A very striking outcome of Fig. 0] 
is that, for integrated quantities as the ones involved in neutrino experiments, the RMSGA 
cross sections compare very well with the RDWIA ones down to remarkably low ejectile 
kinetic energies of about 200 MeV. Below this energy, the RMSGA predictions are not 
realistic due to the underlying approximations, mainly the postulation of linear trajectories 
and frozen spectator nucleons. 

For the sake of completeness, in Fig. |3] we show our predictions for CC //-nucleus cross 
sections. The effects of FSI are of the same order as for NC, and very similar results are 
also obtained within RMSGA and RDWIA down to very low ejectile kinetic energies. 

C. Estimating the effect of FSI mechanisms 

A quantity routinely used to estimate the overall effect of FSI in nucleon-emission pro- 
cesses is the nuclear transparency. Intuitively, it provides a measure for the probability that 
a nucleon of a certain energy - above the particle-emission threshold - can escape from the 
nucleus without being subject to any further interaction. From this 'definition', one can 
expect that the nuclear transparency is identical for neutrino and electron induced nucleon 
knockout. Once the nucleon is traversing the nuclear medium, only its energy is expected 
to determine the way it propagates. In addition, neutrinos and electrons can be expected 
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FIG. 4: Neutral current 12 C(u, v') (left panels) and 56 Fe(i^, u') (right panels) cross sections as a 
function of T/v at different incoming neutrino energies. The solid lines represent the RPWIA 
predictions of the Madrid group, in agreement with those of the Ghent one. The dashed (dot- 
dashed) lines implement the effect of FSI within the RMSGA (RDWIA) framework. 

to probe equal amounts of bulk and surface parts of the target nucleus. 

Several investigations of the nuclear transparency have been carried out using the A(e, e'p) 
reaction in the QE regime (i.e. the Bjorken variable x = Q 2 /{2M^uj) ~ 1), and data for 
different nuclei are now available. The nuclear transparency is extracted from the ratio of the 
measured A(e, e'p) yield to the calculated one using the plane- wave impulse approximation, 
according to 

rp (Q 2 ) — Jy j PmdE m Y exp (E m , p m , kf) 
exp c(A) J v d 3 p m dE m Y PWIA (E m , p m ) ' 

The quantity V specifies the experimental phase-space in missing momentum (p m ) and en- 
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FIG. 5: Charged current 12 C(V M , (left panels) and 56 Fe(f A1 , (right panels) cross sections as 
a function of the outgoing lepton energy e' at different incoming energies. The labeling is the same 
as in Fig. 

ergy (E m ). The kinematics cuts \p m \ < 300 MeV/c and E m < 80 MeV, in combination with 
the requirement that x ~ 1, guarantee that the electro-induced proton-emission process is 
predominantly quasi-elastic. The factor c(A) is introduced to correct in a phenomenological 
way for short-range mechanisms, and is assumed to be moderately target-mass dependent 
(c = 0.9 for 12 C, and c = 0.82 for 56 Fe). It accounts for the fact that short-range correlations 
move a fraction of the single-particle strength to higher missing energies and momenta and, 
hence, beyond the ranges covered in the integrations of Eq. (J52j) . Without going into details, 
theoretical predictions are obtained in a similar way from the ratio of distorted-wave calcu- 
lations to plane-wave ones. In Fig. ® the transparencies predicted by the RMSGA and the 
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RDWIA models are displayed as a function of Q 2 for 12 C and 56 Fe, together with the world 
data extracted from A(e,e'p). Solid (dot-dashed) lines show the A(e, e'p) results within 
RMSGA (RDWIA). Details about the calculations can be found in Ref. 29(. The dashed 
(RMSGA) and dotted (RDWIA) curves correspond to the computed A(v, v'p) transparen- 
cies, obtained using the same procedure as for electron scattering. This procedure includes 
the computation of RDWIA and RPWIA cross sections at x ~ 1, averaged over the same 
phase space used in Eq. (I52j) . As can be seen, within each model the neutrino transparencies 
agree quite well with their electron counterparts. This result clearly illustrates the fact that 
in our models the average attenuation effect of the nuclear medium on the emerging nucleon 
is rather independent of the nature of the leptonic probe. 

Adopting the idea that the nuclear transparency for electrons equals the one for neutrinos, 
the information obtained about nucleon propagation via A(e, e'p) can be used to predict the 
effects of FSI mechanisms in inclusive QE i/-nucleus cross sections. As the transparency is 
essentially the ratio of cross sections including FSI to the ones in the plane- wave limit, this 
will be done by multiplying the RPWIA results for neutrino-nucleus cross sections with the 
measured transparency factors extracted from A(e,e'p). In this scenario, the benchmark 
RPWIA neutrino-nucleus cross sections are crucial. It is important to realize that we use 
transparency factors that are confined to i ~ 1, while the computation of the inclusive 
neutrino-nucleus cross section include the full phase-space. 

In Fig. the dashed and dot-dashed lines represent the inclusive CC ^-nucleus cross sec- 
tion within RMSGA and RDWIA, respectively. The solid curve displays our corresponding 
RPWIA calculation, scaled with a constant factor taken as a representative value for the ex- 
perimental A(e, e'p) transparency for the nucleus. For 12 C ( 56 Fe) we take T « 0.52 (« 0.34). 
In extracting these values, we have corrected the measured transparencies from Fig. |U] with 
the factor c(A). A very good agreement is observed between the rescaled RPWIA and the 
full RDWIA/RMSGA curves in the case of 12 C. This finding supports the idea that a simple 
scaling of the RPWIA results with a transparency factor obtained from electron scattering 
data allows one to reliably estimate the FSI effects for the quasi-elastic contribution to the 
inclusive neutrino cross section. The fact that for 56 Fe the agreement is less satisfactory 
reflects the fact that our models slightly underestimate the 56 Fe transparency data. 

Finally, Fig. |H1 displays the total cross section (er = J de' (da j de')) for 12 C(^ M , ji~) and 
56 Fe(z/ At , fi~) reactions, scaled with the number of neutrons in the target. Results are shown 
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FIG. 6: Nuclear transparencies versus Q 2 for different nuclei in quasi-elastic kinematics. The solid 
(dot-dashed) lines shows the results of a RMSGA (RDWIA) A(e, e'p) calculation 3|- Tne dashed 
(dotted) lines represent the results for A{v,v'p) within RMSGA (RDWIA). Data points are from 
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49| (solid triangles), and 



(diamonds) . 



within RPWIA and RDWIA using a complex optical potential. The figure clearly shows that 
the difference between RPWIA and RDWIA cross sections is approximately given by the 
experimental transparency factor extracted from A(e, e'p) at QE kinematics. Furthermore, 
other important features can be extracted from this figure. First, the RPWIA cross sections 
scale with the target mass-number. In this way, when RPWIA cross sections are required 
for a heavy nucleus, a very good approximation consists in multiplying this cross section 
per nucleon by its mass number. Second, the cross sections do not appreciably change from 
neutrino energies above 2 GeV, i.e. the cross sections saturate at high incoming neutrino 
energies. To finish with, we compare our relativistic calculations with data from various 
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FIG. 7: Charged current 12 C(z^, (left panels) and 56 Fe(^ M , (right panels) cross sections as a 
function of e' at different incoming energies. The dashed (dot-dashed) lines represent the RMSGA 
(RDWIA) prediction. The solid lines show the RPWIA results, scaled with a transparency factor 
T( 12 C) « 0.52 and T( 56 Fe) « 0.34. 

experiments. The RPWIA calculations give a fair account of the neutrino-energy and mag- 
nitude of the data. The RPWIA is confined to single-nucleon knockout thereby not including 
final-state interaction effects. The RDWIA calculations, on the other hand, including FSI 
effects via the introduction of an optical potential, considerably underestimate the data. 
The results contained in Fig. |H] indicate that at least 50% of the measured (i^, fi~) strength 
can be attributed to single-step ("elastic") nucleon knockout to missing energies below 80 
MeV in the residual A-l nucleus. The remaining fraction of about 50% could be attributed 
to multi- nucleon knockout, pion production with subsequent reabsorption, single-nucleon 
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FIG. 8: Total CC (z^,^ - ) neutrino cross sections as a function of the incoming neutrino energy. 
The solid (dashed) line shows the RPWIA calculations on 12 C ( 56 Fe). The dot-dashed (long-dotted) 
curves implement the effect of FSI on 12 C ( 56 Fe) within RDWIA. The short-dotted (long-dashed) 



show the RMF curve for C ( Fe). All results are sea 



target. Data points are from Refs. 
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Adding all these contributions would move the 
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431 ] that the importance of these missing channels can be 



estimated in a model in which both the single-particle bound-states and scattering states 
are computed in a real mean-field potential obtained in the Hartree approximation to the 
<t — uj model. Such an extreme mean-field model (here coined as RMF) involves no imagi- 
nary potential and accordingly the loss of single-nucleon knockout strength into "inelastic" 
nucleon-knockout channels is effectively reintroduced. As can be appreciated from Fig. |Hl 
the presence of a real potential reduces the RPWIA single-nucleon knockout strength, but 
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the global reduction with respect to RPWIA is much smaller compared to the 50 — 60% 
in the RDWIA/RMSGA frameworks. It is remarkable that the energy dependence of the 
cross section is identical in all three relativistic frameworks adopted here. The magnitude, 
on the other hand, depends strongly on the model used to account for FSI mechanisms. 
This reinforces our suggestion that one could use RPWIA results to predict the "elastic" 
single-nucleon knockout contribution to inclusive neutrino cross sections, provided that they 
are rescaled with a transparency factor extracted from A(e, e'p) data. 

IV. CONCLUSIONS 

We have employed two relativistic models, RMSGA and RDWIA, to study NC and CC 
quasi-elastic neutrino-nucleus scattering. Results have been presented for carbon and iron 
targets, covering a wide range of neutrino energies. Within RPWIA both models provide 
nearly identical results, which deviate from existing RPWIA predictions. The fact that 
two independently developed codes that adopt very different numerical techniques agree in 
this limit, together with the fact that our RPWIA predictions approach the RFG model 
at high energies, give us confidence that our RPWIA calculations serve as benchmark re- 
sults. We subsequently computed the effects of FSI mechanisms within the RMSGA and 
RDWIA models. The two ways of dealing with FSI are consistent down to remarkably low 
outgoing nucleon kinetic energies of about 200 MeV. FSI produce a large reduction of the 
cross sections, that increases with the mass number of the target nucleus. Finally, we have 
illustrated that the nuclear transparencies extracted from A(e, e'p) measurements can be 
used to estimate the effect of FSI mechanisms on the elastic single-nucleon knockout con- 
tribution to the inclusive neutrino-nucleus cross sections. Extensions of our models include 
the implementation of pion production and delta resonance. Work in this direction is in 
progress. 
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